Paramagnetic properties of non-centrosymmetric superconductors: 

Application to CePtaSi 
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In a non-centrosymmetric crystal, the Zeeman interaction of the band electrons with an external 
magnetic field is highly anisotropic in the momentum space, vanishing along some high-symmetry 
planes. One of the consequences is that the paramagnetic susceptibility in superconductors without 
inversion symmetry, such as CePtaSi, shows an unusual temperature dependence. 
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Electronic spin susceptibility measurements in the su- 
perconducting state provide one of the most useful tools 
for the identification of the pairing symmetry. Exper- 
imentally, the susceptibility x(T) is measured by the 
Knight shift in the nuclear magnetic resonance frequency 
due to the hyperfine interaction of the conduction elec- 
trons with the nuclear magnetic moments. If the su- 
perconducting pairing occurs in the spin-singlet channel, 
then the Cooper pairs have spin S = Q and do not con- 
tribute to the magnetization M of the system, which is 
therefore entirely determined by the thermally excited 
quasiparticles. For a fully gapped order parameter one 
would see an exponentially decreasing x{T) at low tem- 
peratures. In contrast, if the gap has zeros at the Fermi 
surface, then x(r) oc for isolated point nodes, or 
x(T) oc T for line nodes. In the triplet case, when the 
Cooper pairs have spin 5=1 and the order parameter 
is a vector d{k) in the spin space (see, e.g., Ref. 0), the 
susceptibility depends on the mutual orientation of d and 
the external magnetic field B. If B || d, then the Cooper 
pairs do not contribute to M , and x(T) has the same 
temperature dependence as in the singlet case. If B _L d, 
then both the pairs and the excitations contribute to M , 
so that x(^) = Xn ~ the normal-state susceptibility. The 
observation of a flat Knight shift in Sr2Ru04 for the field 
in the basal xy plane 2], has been used as a proof of a 
spin-triplet pairing with d || z. 

The theoretical picture described above is valid only 
if the superconducting crystal has an inversion center. 
Although this is the case in the majority of superconduc- 
tors, some exceptions have been known since 1960's 3]. 
More recently, it was pointed out in Ref. |1| that the 
surface superconductors, e.g. Na-doped WO3 [3, are in- 
trinsically non-centrosymmetric simply because the two 
sides of the surface layer are manifestly non-equivalent. 
As for the bulk materials, the latest examples are CePtsSi 
@ and Ulr 0. Different models of superconductivity 
in CePtaSi have been proposed in Refs. [1,0. In this 
Letter, we calculate the suppression of the critical tem- 
perature and the paramagnetic susceptibility in super- 
conductors without an inversion center. We focus on the 
tetragonal symmetry relevant for CePtaSi, and show in 
particular that the anisotropy of the temperature depen- 
dence of the susceptibility tensor Xij (T) is strikingly dif- 



ferent from the centrosymmetric case. 

The spin susceptibility in two-dimensional non- 
centrosymmetric superconductors has been previously 
studied in Refs. |3. lift llll ll ^ . where the inversion sym- 
metry breaking in the presence of a non-zero spin-orbit 
(SO) coupling was introduced using the Rashba model 
[l3| . Very recently, a three-dimensional generalization of 
the Rashba model was applied to CePtsSi in Ref. [l^ . 
Treating the SO band splitting as a perturbation, it was 
found that the order parameter becomes a mixture of 
a spin-singlet (even in k) and a spin-triplet (odd in k) 
components, which gives rise to a non-zero residual sus- 
ceptibility at T = 0. In this Letter, we use a differ- 
ent approach based on the effective single-band Hamil- 
tonian, which works for any crystal symmetry and arbi- 
trary strength of the SO coupling. Guided by the fact 
that the SO band splitting is usually large compared to 
the superconducting energy scales, we consider both the 
Cooper pairing and the magnetic response independently 
in different bands. In contrast to the previous works, we 
construct the pairing interaction using the exact band 
states and explicitly take into account that all the pair- 
ing channels but one are suppressed [l5| . 

The starting point of our analysis is the observation 
that in a crystal lacking an inversion center the elec- 
tron bands are non-degenerate almost everywhere, except 
along some high-symmetry lines in the Brillouin zone. 
Indeed, without the inversion operation /, one cannot 
construct two orthogonal degenerate Bloch states at the 
same k. At zero SO coupling there is an additional sym- 
metry in the system - the invariance with respect to ar- 
bitrary spin rotations, which restores the two- fold spin 
degeneracy of the bands. Here we assume that the SO 
coupling is sufficiently strong, so that the bands are well 
spht. The results of Ref. "^I show that this is indeed the 
case in CePtsSi, where the SO band splitting can be as 
large as 50-200 meV depending on the band. Assuming 
that there is no disorder in the crystal, the Bloch wave 
vector fc is a good quantum number in zero field. The 
free electron Hamiltonian for a non-degenerate band can 
be written as Hq = J2k^i^)'^k'^k, where J^k stands for 
the integration over the first Brillouin zone, and e(k) is 
the quasiparticle dispersion, which takes into account the 
periodic lattice potential and the SO interaction. 
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If the time-reversal symmetry is not broken in the nor- 
mal phase, then the states |fc) and K\k) ~ | — fe) are 
degenerate because of the Kramers theorem 16|. It is 
the coupling between those states that leads to the for- 
mation of the Cooper pairs and the superconductivity in 
the system. The large band splitting strongly suppresses 
the pairing of electrons from different bands. Then, con- 
sidering just one band (i.e. neglecting the inter-band pair 
scattering) and assuming that the interaction has a gen- 
eralized Bardeen-Cooper-Schrieffer (BCS) form, we have 

H„u ^\Y. V{k,k')clcl^c^k'Ck'. (1) 

k,k' 

The pairing potential can be written as 
Vik.k') = V{k,k')t{k)t*{k'), where V{k,k') ^ 
—Vr ^a{k)(j)*a{k') is the part that transforms accord- 
ing to an irreducible representation T of the normal-state 
point group G, (l)a{k) are the scalar basis functions of 
r, which are are nonzero only inside the energy shell 
of width Ec near the Fermi surface, Vr > is the 
coupling constant, and t(k) = —t{—k) are non-trivial 
phase factors in K\k) = t{k)\ — k) Although 
anti-commutation of fermionic operators dictates that 
the mean-field order parameter A(fe) — t{k) ria4>a{k) 
is odd in fc Q, its nodal structure is determined by the 
basis functions 0a (fc), which should be even because of 
the presence of t{k) jlSj. The focus of this article is on 
CePtaSi, which has a non-centrosymmetric tetragonal 
crystal lattice described by the point group G — C^y. 
This group is generated by the rotations about the 
z axis by an angle 7r/2 and the reflections ax in the 
vertical plane (100), and has five irreducible represen- 
tations: four one-dimensional (ID): Ai, A2, Bi, B2, 
and one two-dimensional E [T9l |. Here are the examples 
of the even basis functions: oc ~^ ~^ ^/c^, 

4>A2 oc kxky(kl ~ ky), (pBt cc kl ~ ky, 4>b2 oc k^ky, 
and {4>E,i,4'E,2) oc [k^kz, kyk^). We consider a small 
sample, of a dimension d < ^ < S, where ^ is the 
superconducting correlation length and S is the London 
penetration depth, which allows us to neglect the spatial 
variations of both the order parameter components rja 
and the magnetic field. 

Let us now turn on a uniform stationary magnetic 
field B = curl A. Assuming that the pairing interac- 
tion is field- independent, B can only affect the system 
through its couphng to the band states. At B ^ 0, the 
band dispersion function e{k) is replaced by an effective 
band Hamiltonian in the momentum space, which can be 
represented as a power series in B: e(fc) — » £{k,B) = 
e{K) + Bieis{K) + where K = k + {e/hc)A(iVk) be- 
cause of the requirements of gauge invariance |2(]| | . The 
expansion coefficients must satisfy certain symmetry- 
imposed conditions, in particular the zero-field band dis- 
persion e(fc) must be invariant under all operations from 
G. In addition, at B 7^ the Hamiltonian is invariant 
with respect to time reversal K only if the sign of B 
(and of A) is also changed, which imposes the following 



constraint on the function S: K^£{—B)K = £{B). 

In the analysis of the "paramagnetic" properties of su- 
perconductors, the orbital effect of the field is neglected, 
which is achieved by putting A = in the effective band 
Hamiltonian. Then, for a two-fold degenerate band in the 
presence of inversion symmetry, f is a 2 x 2 matrix, and 
the coupling to the magnetic field is described by a famil- 
iar Zeeman term: £ap{k,B) = e{k)5ai3 — Bi^ij {k)a j^a/s, 
with idij{k) = iiij{—k) being the tensor generalization of 
the Bohr magneton /xb for the case of band electrons. 
The indices a (3 here are pseudospin indices ,,21j . The 
Zeeman interaction splits the energies of the electrons 
forming the Cooper pairs and gives rise to the paramag- 
netic suppression of superconductivity [2^ . 

If the inversion symmetry is absent and the bands are 
non-degenerate, then the Zeeman term should be modi- 
fied. The effective single-band Hamiltonian in the exter- 
nal field can be written as 

Ho = J2Hk)-BXik)]cick, (2) 

k 

which is markedly different from the centrosymmetric 
case. Here \{k) is a real pseudovector, which satisfies 
the conditions {gX){g~'^k) = A(fc), where g is any op- 
eration from the point group G. Because of the time- 
reversal symmetry, we also have e(— A;) — e{k) and 
A(-fc) = -A(fc), but £{-k, B) ^ £{k, B). 

Explicit expressions for A(A;) can only be obtained 
in some simple models. For example, in an isotropic 
two-dimensional electron gas in the xy plane with G = 
Coon, the combined effect of the SO coupling and the 
lack of inversion symmetry is described by an addi- 
tional (Rashba) term in the single-particle Hamiltonian: 

Hso = lJ2k'"' ■ ('^'T't' ^ k)al^aka', where n is the 
normal vector to the plane Diagonalization of 

the Hamiltonian in zero field gives two non-degenerate 
Rashba bands: e±{k) = eo{k) ± 7|fc|. At finite field, 
adding a usual Zeeman term ~^b{B ■ cr), and expand- 
ing the eigenlavues of the Hamiltonian in powers of B, 
we obtain £±{k, B) = e±(fe) - \±{k)B + 0{B^), where 
A±(fe) = ±fiBik X n)/\k\. Thus the coupling of the 
Rashba bands with the field is highly anisotropic, in par- 
ticular it vanishes for S || n. 

While a microscopic derivation of the effective single- 
band Hamiltonian Q in more realistic systems can be 
done, at least in principle, using the procedures described 
in Refs. |23|, it suffices for our purposes to work with 
a phenomenological expression for \{k), which is com- 
patible with all the symmetry constraints. We need an 
expression for X{k), which satisfies (i) A(— fc) = — A(fe), 
(ii) (C4.A)(C4-/fc) = A(fc), and (iii) (a,A)((7-ifc) = A(fc) 
(since A is a pseudovector, we have ct^jA = IC2xX = 
C2xA, where C2X is a rotation by an angle tt about the 
X axis). It is straightforward to check that the general 
expression for A(A;) is given by 

A(fc) = (j)E.2{k)S: - (t)E,i{k)y + 0^2(^)2, (3) 
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where (f>E,i{2){f^) E^nd (/^^^(fc) are real odd functions, 
which transform according to the representations E and 
A2 respectively, e.g. (f>A2 oc k^kykzik"^ — ky) and 

{4>E,i,4>E,2) oc {kx, ky). We see that Xz = along the five 
nodal planes of the A2 representation, while A^: = Aj, = 
along the z axis. 

Now, we calculate the free energy T for the Hamilto- 
nian H = Ho + H^nt, defined by Eqs. ^ and l^l. We use 
the effective field theory in terms of the bosonic Matsub- 
ara fields »7a(T), which can be introduced in a standard 
fashion by decoupling the pairing interaction The 
effective action for a uniform stationary order parameter 
in the mean-field approximation reads 



where 



- e{k) + B\{k) 



ixrlnG-i, 



-A(fc) 

e{k) + BX{k) 



(4) 



(5) 



is the 2x2 inverse Gor'kov Green's function, w„ = 
(2n + l)7rT is the fermionic Matsubara frequency, and 
V is the system volume. The mean- field free energy is re- 
lated to the saddle-point action T — (l//3)S'e//, and 
the mag netization density is M = -V^'^idT/dB). The 
saddle point condition yields the self-consistency equa- 
tion 



(6) 



which determines the temperature and field dependence 
of the order parameter components. Substituting here 
the Green's function |SJ, one can see that the phase fac- 
tors t{k) drop out of the gap equation. 

Let us first find how the critical temperature is sup- 
pressed by the field. The equation for Tc{B) is obtained 
by linearizing Eq. JH)) and can be written in the form 
det llifabll = 0, where 



In- 



't'c 



(7) 



0:(fc)0,(fc)ReV 



.-BA(fc) 

' 27rT, 



FS 



where ^{x) is the digamma function, T^o — 
l.lScc exp(l/iV^FT^) is the critical temperature in zero 
field {Np is the density of states at the Fermi level), and 
the angular brackets denote the average over the Fermi 
surface. Next, differentiating the effective action with 
respect to B, we find the magnetization density 



M = ^ ^ A(fc) [Gn (fc, a;„) + G22 (fe, ojn) 



(8) 



n k 



The uniform susceptibility tensor is defined in the usual 
manner as Xij = dMi/dBj\B=Q- One can easily check 



using Eq. © that the corrections to the order parameter 
components rja in a weak magnetic field are quadratic in 
B, which means that in the calculation of Xij one can 
neglect the field dependence of A(A;), to obtain 



A.(fc)A,(fc) 
4T ^ cosh^[£;(fc)/2T] 



(9) 



where E{k) = y/e'^ik) + |A(fc)|2. 

We now apply the general theory to CePtaSi. The 
pseudovector A(fe) in this case is given by Eq. Q. Since 
all three components of A are in general non-zero, we 
expect the superconducting critical temperature to be 
suppressed for all orientations of the magnetic field. Ac- 
cording to Eq. (0, the magnitude of the suppression 
depends on many factors: the shape of the Fermi sur- 
face, the symmetry of the order parameter, and also the 
explicit form of the functions (j> in Eq. (PJ- The Fermi 
surface of CePtaSi is quite complicated and consists of six 
sheets Q . It is not known which one (or ones) of them 
are superconducting. The order parameter symmetry is 
not known either, although the observation of a linear 
T-dependence of the specific heat @ probably indicates 
that the order parameter has lines of nodes. In view of 
all this uncertainty, it seems to be premature to discuss 
the paramagnetic suppression in CcPtaSi quantitatively. 

More interesting qualitative conclusions can be drawn 
from the analysis of the susceptibility Xij {T) ■ It follows 
from Eq. Q that the main contribution to the suscep- 
tibility tensor at low temperatures comes from the ther- 
mally excited nodal quasiparticles. If the order parame- 
ter A(fc) has no zeros at the Fermi surface (e.g. for the 
Ai representation), the susceptibilities are exponentially 
small in all directions. Let us consider the ID order pa- 
rameters corresponding to the representations A2, Bi, or 
B2, for which the lines of nodes are symmetry-imposed. 
In this case, the components of the susceptibility tensor 
are given by Xxx = Xyy = X\\ and Xzz = X±, where 



XII m 



X±(T) 



-E 

ST ^ 



1,1 ( 



^1,2 (fc) 



8T^ cosh'[£:(fc)/2r] 



1 

4T 



E 



cosh^[E{k)/2T] 



(10) 
(11) 



It is straightforward to show that for the field in the basal 
plane, X\\{T) oc T, since 0|; ^{k) + ^i^) is in general 
non-zero everywhere, except the poles of the Fermi sur- 
face. This behavior is characteristic of the systems with 
lines of nodes. In contrast, for the field orientation along 
the z axis, we have from Eq. Ullfl x^{T) , since 

(t>A2{k) vanishes at the nodal lines for all the ID order 
parameters. Such temperature dependence is never seen 
in the centrosymmetric case. For a two-component order 
parameter A(fc) ~ rjik^kz +r]2kykz with a horizontal line 
of zeros at kz = 0, one also has a T^-behaviour of xx(T). 
On the other hand, the anisotropy and the temperature 
dependence of XnC^) is determined by (?7i,?72). 
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We would like to note that, in contrast to Refs. P.[T^. 
our approach does not yield a finite value of the suscepti- 
bility at T = 0. The explanation is that the residual sus- 
ceptibility in the systems described by the Rashba model 
comes from the interband transitions between the Rashba 
bands (Van Vleck susceptibility), which are not affected 
by the transition into the superconducting state 'l2!|. In 
our model, the quasiparticles in each band respond to the 
external field independently of the other bands. If to take 
the interband transitions into account, then the observed 
susceptibility will be Xtot{T) = Xo + x(2^); where the first 
term is the temperature-independent background that 
comes, e.g., from the Van Vleck processes, and the second 
term is the single- band contribution jnj. Other mecha- 
nisms that might affect the residual susceptibility include 
the contributions from the unpaired sheets of the Fermi 



surface, or spin-reversing scattering at impurities or sur- 
face imperfections p^ . 

In conclusion, we have shown that the paramagnetic 
responses of superconductors with and without inversion 
center are qualitatively different. The most important 
feature is that, in the latter case, the coupling of the 
non-degenerate band electrons with the external field is 
strongly momentum-dependent and vanishes, for sym- 
metry reasons, along some high-symmetry planes in the 
Brillouin zone. This results in a high anisotropy of the 
susceptibility in the superconducting state with lines of 
nodes, from X\\{T) oc T to X-l(^) °c T^, which can be 
used as a clear-cut experimental test of our theory. 

This work was supported by the Natural Sciences and 
Engineering Research Council of Canada. 
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